A Renormalization Group Analysis of Coupled Superconducting and Stripe Order in 

1+1 Dimensions 
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In this paper we perform a renormalization group analysis on the 1+1 dimensional version of a 
previously proposed effective field theory |J describing (quantum) fluctuating stripe and supercon- 
ductor orders. We find four possible phases corresponding to stripe order/disorder combined with 
superconducting order/disorder. 
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I. INTRODUCTION 

In La2- x Sr x CuOi (LSCO) compounds, there are 
three well-established ordering tendencies: antiferromag- 
netism, superconductivity and charge/spin stripes. JjJ. 
Some experiments indicate that stripes and superconduc- 
tivity can even coexist in these compounds P[. Further- 
more, neutron scattering experiments by Lake et al.^ 
show that a moderate magnetic field can have large ef- 
fects on the incommensurate magnetic fluctuations. This 
is widely taken as evidence suggesting the stabilization 
of stripes by the magnetic field. 

In mean field theory, when two order parameters are 
in close competition it is possible for them to coexist in a 
certain region of the phase diagram[Q. In such a coexis- 
tence region quantum fluctuations of both order param- 
eters dominate the low-energy physics. In a recent paper 
LccQ examined such a situation. That paper described 
how the Goldstone modes of stripe and superconducting 
orders and their respective topological defects interact. 

We stress that the theory presented in Ref . Q differs in 
important ways from the conventional self-dual charge- 
density wave/superconductivity action in one dimension. 
Indeed, in one dimension the displacement field of the 
charge density wave is conjugate to the phase of the su- 
perconducting order. As a result the charge density wave 
and the superconducting orders are mutually exclusive 
(i.e. whenever superconducting susceptibility strongly 
diverges the charge density wave susceptibility does not 
and vice versa). In contrast, in the theory of Ref. || there 
exists a generic region in phase diagram where both or- 
ders exist. 

In this paper we examine in detail a one-dimensional 
analog of the model studied in Ref. Q . The motivation for 
this is that well-developed calculational methods (such 
as the renormalization group) can be used to analyze the 
phase structure of the model. This can be used to check 
the correctness of the asserted phase structure in Ref. || 

Now we describe the theory proposed by Lee [|| . Since 
the stripe order is a one dimensional charge density wave, 
its Goldstone mode (i.e. stripe displacement) is a U(l) 
scalar]^. The superconducting order, of course, also pos- 
sesses a U(l) Goldstone mode. The important question 
is: how do these two U(l) modes couple together? A hint 



of how this coupling works comes from the experimental 
fact that the period of incommensurate spin correlation 
decreases as the doping density increases. Motivated by 
this LecH] constructed the following Lagrangian density 
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In the above <j)Q — e l9s is the U(l) phase factor of 
the superconducting order parameter, and po — e l6p 
is the phase factor of the stripe order. That is, 9 p = 
(2tt/A)x- u(x, t), with A the stripe period and u(x,t) the 
displacement field of the stripe order. J M and are aux- 
iliary fields coupling to the superconducting and stripe 
phases, respectively. These auxiliary fields have the phys- 
ical interpretation of energy-momentum currents. (In 
this paper greek indices run over x, t and repeated in- 
dices are summed.) 

Without the coupling {gx,i = 0), integrating out J p 
and q^ produces the field theory for two independent 
U(l) Goldstone modes and their respective vortices. The 
effect of the coupling is to favor stripe displacement (u) 
in the presence of local charge imbalance (Jt). 

To analyze Eq. (|l]), Lee used a duality transformation 
plus an educated guess about the four possible quantum 
phases corresponding to combinations of stripe and su- 
perconducting order/disorder. In this paper we study 
a one-dimensional version of Eq. (^), applying the well- 
developed techniques of duality transformation and the 
renormalization group to determine the possible phases 
in a more unbiased fashion. We find that all four com- 
binations of stripe order/disorder and superconducting 
order/disorder are stable phases. This supports the con- 
jectured phase structure in Ref. I|. 

In the following we use a real-space renormalization 
procedure similar to that used by Kosterlitz and Thou- 
less to treat the phase transition of the 2-dimensional 
coulomb gas@, ||. In section two we derive the vortex 
action (the vortex of the stripe order parameter is the 
dislocation). In section three we obtain the renormaliza- 
tion group recursion relations for the coupling constants 
in that theory. As in the Kosterlitz-Thouless theory we 
make the small vortex fugacity approximation. We ana- 
lyze the implications of these flows for phase stability in 
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section four. 



II. DUALITY TRANSFORMATION TO 
2-SPECIES COULOMB GAS 

Following the work of Jose et al|| we first perform a 
duality transformation and write the theory in terms of 
vortex degrees of freedom. 

Starting with Eq. (|l|), we first separate the phase of cf>o 
and po into a topologically trivial part and a topologically 
non-trivial part: 



Po = e*°p 



(2) 



In the above r]o and £o are single valued, while <f> and 
p contain configurations with non-zero windings. After 
integrating over the topologically trivial phases (?7o,£o) 
we obtain two conservation laws: 



= 
= 0. 



(3) 



To explicitly fulfill these conservation laws we write J M = 
t^ v d v A and = e^d^x, where x an d A are scalar fields. 
Substitution leads to 
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~ , )T ~ -(d^xf + T^id^f + e^dMd^ 

+e^d u xpdf,p + igidtxdx^ + igAxdtA. (4) 

Upon integrating by parts and identifying the vortex den- 
sities N — ie-^udvipd^p) and M — ie^d^^d^) the La- 
grangian density becomes 
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(d^A) 2 + iAM + i X N 



(5) 



The above equation can be written in momentum space 
as 



+ i(x(k)A(k))" 



(N(k)\ 
{M(k)J 



(6) 



where G — g\ + g 2 . Integrating out the x an d A fields 
then produces 
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l(NM)* — 
2 V 'det 
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(8) 



Eq. (||) is the starting point of our renormalization 
group analysis. It describes a system of two interacting 
(anisotropic) coulomb gases - the vortices of the super- 
conducting order parameter and the dislocations of the 
stripe order parameter. Inspired by the work of Koster- 
litz and Thoulcss we perform a real space renormalization 
group analysis of Eq. m) in the following. 



III. RENORMALIZATION GROUP ANALYSIS 

Eq. (||) is more complicated than the one species 
Coulomb gas problem in two respects: 1) there are two 
species of vortices and 2) the interactions are not rota- 
tionally invariant (i.e. the interaction depends not only 
on the distance between vortices but also on their relative 
orientation). In order to complete the renormalization 
group program we have to characterize the interaction 
in terms of a discrete set of coupling constants. One 
way of achieving this is to Fourier transform the angular 
dependence of the vortex-vortex interaction. In momen- 
tum space each element of the interaction matrix is of 
the form G(k) = G(fc, 9) = sMl (Q is the angle made by 
k and the k x axis) . Therefore we expand each of these 
terms in a fourier series, e.g. g{9) = J2 n a ne' m9 '- When 
transformed back to real space, our action then becomes 



S = ~ y*d 2 R 1 d 2 R 2 iV(R 1 )G 7V (R 1 -R 2 ) N(R 2 ) 



+Af(R 1 )G M (R 1 -R 2 )M(R 2 ) 
+2M(R0 ir(Ri - R 2 ) iV(R 2 ), 



where 



Gjv(R) 



G M (R) = 



21, / ~ 



r(R) 



In the above 
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(9) 



(10) 



(11) 



(12) 



\n— — oo 



a n = (-l) n a*_ n 
a n = (-1J a_ n 



(-i)"ci„ 



(13) 



to ensure the interaction functions are real. We stress 
that because of the angular dependence of Eqs. [n| O , 
and [l^, Gn,Gm and T depend not only on the dis- 
tance |Ri — R2I but also on the relative orientation 
(Ri — R 2 )/|Ri — R2I. In general a n and a n are nonzero 
only for even n. The physical reason for this is indistin- 
guishability of two charges of the same type (for details 
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see the appendix). c„ can be nonzero for both odd and 
even n. 

The limit where all the a n , a n , and c„ are zero except 
ao and ao describes two decoupled isotropic 2-D X-Y 
models in their coulomb gas representations - the Koster- 
litz and Thouless problem. Before we attack Eq. (@), as a 
warm up, let us briefly review the Kosterlitz-Thouless re- 
sults for the one component system. In the renormaliza- 
tion group approach one integrates out one pair of tightly 
bound dipole (i.e. a dipole with r c + dr c < size < r c ) at a 
time. The renormalization group proceeds iteratively by 
treating r c as a running length scale. The two coupling 
constants in that case are the vortex-vortex interaction 
strength ao and vortex fugacity y — e _Al where \x is the 
core energy of vortices. In the limit of y << 1, the renor- 
malization group equations for ao and y are given by 



dy 

dl 
dag 

~dT 



(>-£) 



2 2 



(14) 
(15) 



The above equations have the entire y = axis as fixed 
points. However depending on whether ao — 8ir is pos- 
itive/negative the fixed point is stable/unstable. The 
point y = 0, ao = 87r is a critical point. Near it, the flow 
trajectories are given by 



(4tt)V 



C 



(16) 



Here C is a constant labelling each trajectory. This flow 
is shown in figure [|. Note that the C = separatrix y = 
(4tt) 4 ( a ° — ^ 7r ) separates the basins of attraction for the 
ordered and disordered phases. In the ordered phase the 
density of vortices renormalizes to zero (y — > 0) at large 
length scales, signifying the presence of a bound dipole 
phase. In the disordered phase the density of vortices 
increases (y increases) at large length scales, signifying 
the existence of a vortex plasma (unbound dipole) phase. 

The two-component vortex gas problem we are fac- 
ing is not so different. However, we have to painstak- 
ingly keep track of all the Fourier coefficients in Eq. (10), 
Eq. dl) and Eq. and examine how they renormalize. 
Interestingly, even in the presence of these anisotropic in- 
teractions the Kosterlitz-Thouless renormalization group 
program closes. 

All the technical details are given in the appendix. 
Here we just note the following point. Because the 
"Coulomb charge" of the superconducting vortices is not 
related to the "Coulomb charge" of the stripe vortices, 
only intra-species dipoles are possible. This implies that 
the positions of vortices belonging to different species do 
not have to obey the constraint that the minimum dis- 
tance is r c . To lowest order in y the resulting renormal- 
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FIG. 1: Renormalization group flow for the Kosterlitz- 
Thouless transition. The fixed points are the y=0 axis. In 
our model this corresponds to only y, ao nonzero. 



ization group equations are given by 
dyN 



dl 
dyM 
dl 
da n 
~~dT 

da n 
~df 

dc n 

~~dT 



UN 



(»-£) 



VM 2 



Oo 
47T 



-7T ^ S n>k+m (y 2 N a k a m + y 2 M (-l) m+1 c k c m ) 

k,m 

-7T ^ S n ,k+m (VMOtka m + y%(~l) m+1 C k C m ) 
k,m 

n,k-\-m (VN a k + VM<Xk) Cm. (17) 



Note that the renormalization of yjv.M only depends 
on the isotropic part of the coupling (ao and ao). The 
renormalization of an anisotropy coefficient, (e.g. a n ) 
includes many terms. Each term is a quadratic function 
of a*, a/c, or Cfe. If we set all coupling constants except 
ao and ao to zero we recover the Kosterlitz-Thouless flow 
equations (for two separate species). It is easily verified 
that the condition for Gn, Gm, and T to be real (a„ = 
(— l) n al n , etc.) is preserved by these flow equations. It 
is also clear from the form of these equations that these 
coefficients form a closed set under renormalization. 



IV. PHASES OF THE TWO-SPECIES 
COULOMB GAS 



Eq. ( |l7j ) predicts fixed points for yM — yN — and 
a n , a n , c n — anything. As in the normal Kosterlitz- 
Thouless case, we interpret y = as the absence of un- 
bound dipoles. From first line of Eq. (jl7|), y^ = is 
linearly stable if ao > 87r. For ao > 87r the renormaliza- 
tion group brings yjv to larger values. Similar statements 
hold for ao and yM- This suggests the presence of four 



4 



phases depending on whether the vortices of N or M 
species form dipoles or unbind. 

However, this is not quite enough for our purposes. 
What we really need to know is whether all four phases 
can be reached by varying the five parameters in Eq. (|7]) . 
Put another way, the physical system of Eq. (R) lives in 
a five dimensional subspace of the infinite dimensional 
space formed by the a„, a„, and c„. We need to check 
which phases can be reached by trajectories originating 
in the physical subspace, not just which phases exist for 
the infinite dimensional space. 

In order to obtain a tractable problem, in the follow- 
ing we concentrate on the case in which K pp = K p , 
K$ p = K,f, and GK p ^ « 1. When GK P:< p are small 
it is easy to evaluate the fourier coefficients in equations 
Eqs. §,0 and |l|in powers of GK p ^. If GK p< p = 0{e) 
we find that the leading contribution to a„, a„, c„ is 
C(el"l/ 2 ). For the specific form of interaction in Eq. (0) 
it is simple to see that besides Equations |l3| there are 
additional constraints on a n ,a n and c n : 



o-n = &n — unless n = Am 
On = unless n — 4m + 2. 



(18) 



All of these conditions are preserved by the flow equa- 
tions. In terms of the original parameters in Eq. (|7|), the 
non-vanishing coefficients up to order e are: 



a = K p (l- K^K p G 2 /8) 
a Q = K^(l- K^K p G 2 /8) 
c 2 = c*_ 2 = iGKpK^/4 



(19) 



In the following we truncate the space considered to only 
these coefficients, which is correct to lowest order in e. 
Furthermore, this restricts us to a five-dimensional space 
of parameters, which we can take to be independently 
determined by the five parameters in Eq. (|7|). 
In this case, the flow equations [l?] become 



dyN 

dl 
dyu 

dl 

ddQ 

~dT 

dctQ 

~dT 
dc 2 
~df 



VN 



(»-£) 



-K(y 2 N al ~y 2 M \c 2 \ 2 ) 

-n {VM a o -Vn\ c 2\ 2 ) 
-7T (y 2 N a Q c 2 + y 2 M a Q c 2 ) 



(20) 
(21) 
(22) 
(23) 
(24) 



First, by multiplying Eq. (|24| ) by c 2 and then adding 
the result to its complex conjugate, we obtain 



d\c 2 \ 
dl 



-2tt (y N a 



C2 



(25) 



Using this in Eq. (|22j) and Eq. (|23|) then gives us 

da-o 2 I 2 , i i2\ 1 d\c 2 \ 

— = -ny N (a + |c 2 | ) - j^—jj- 



dao 
~df 



C2 



1 d\c 2 \ 
16tt dl 



(26) 

(27) 
(28) 



At this point, we examine closely the region of parame- 
ter space around the critical point by making the change 
of variables 

a = ao — 87r 
a = a a — 8tt 
c = c 2 —c 

In the above c is the fixed point of c 2 . After some algebra, 
and keeping terms to lowest order in a, a, c, 2Mr, and yMi 
the flow equations for a and a are, after some algebra, 



dy 2 N 
dl 
dy 2 



1 

2^ 



ylia 



dl 

da^ 
dl 
drf 
dl 



-2ny 2 N a ((8n) 2 + |c| 2 ) 
-27^0 ((Stt) 2 - |c| 2 ) 



a d\c 2 \ 2 
16n dl 

a d\c 2 \ 2 
16tt dl 



Finally, we can combine these equations to yield 



Xy 



N) 



dl 

d{a 2 ~Xy 2 M ) 
dl 



dC N 

~df z 

dCni 
' dl 



d\c 2 \ 



8tt dl 

a d\c 2 \ 2 
~8^ dl 



(29) 

(30) 

(31) 

(32) 
(33) 

(34) 
(35) 



Here X = (4tt) 4 (1 + |c| 2 /(8tt) 2 ). To understand these 
equations, note that the quantity in parentheses on the 
LHS of each equation is precisely of the form of the con- 
tour numbers for trajectories in the Kosterlitz-Thoulcss 
case (cf. Eq. (|l6|)), with the slope of the separatrix renor- 
malized from -jj^yi to 1/X. In the absence of interac- 
tions (c = 0) these contour numbers Cn and Cm are 
conserved, but in the presence of interactions the renor- 
malization group pushes the flow from one contour to the 
next. Furthermore, by Eq. ( |34| ) and Eq. (|2|), for a > 
(i.e. ao > 8tt) the contour number C'n increases, while 
for a < (i.e. a < 8tt) Cn decreases. Similarly for Cm 
and a. The resulting flow is diagrammed in figure ^[ 

It is clear from this that for a trajectory originat- 
ing below the separatrices (jjn < (ao — 87r)/A and 
yu < (c*o — 8n)/X), the flow leads to both y^ and 
yM zero. In other words, there is a stable phase with 
both types of vortices bound as dipoles, corresponding to 
a stripe ordered/superconducting ordered phase. If the 
trajectory starts with (ao — 8ir) < and (a^— 87r) < 0, the 
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a (« ) 



stripe order 
SC disorder 



stripe disorder 
SC disorder 



stripe order 
SC order 



K. 



stripe disorder 
SC order 



FIG. 2: Renormalization group flow for weakly coupled su- 
perconducting(stripe) vortex gases, with contour numbers C 
labelled. The fixed points are the y=0 axis. With nonzero 
coupling, the flow moves from one contour to the next as in- 
dicated by the arrows. It is clear that even in the presence of 
coupling for each gas there are still two phases corresponding 
to 2/ — > and y increasing. 



FIG. 3: Phase diagram for interacting stripe and Supercon- 
ducting order, for isotropic couplings Kj, and K p and weak 
interaction G. Only the — K p plane is shown. Al- 
though different values of G correspond to different long- 
range interactions, G does not control the existence of 
stripe/superconducting order or disorder. 



flow leads to both j/jv and yu increasing and unbound 
dipoles of both species. Thus there is a stripe disor- 
dered/superconducting disordered phase. Finally, in the 
mixed case, e.g. yjq < (do — 8tt)/X and (arj — 87r) < 0, 
Dn flows to zero but j/m increases. Thus there are phases 
with stripe disorder/superconducting order and stripe or- 
der/superconducting disorder. 

We have seen that, with weak (gi + g%) all four 
phases corresponding to stripe order/superconducting 
order, stripe disorder/superconducting order, stripe 
order/superconducting disorder, and stripe disor- 
der/superconducting disorder are stable and can be re- 
alized in the system described by Lee's theory in (1+1) 
dimensions (Eq. (|l|)). We emphasize that we have ana- 
lyzed only the case with isotropic couplings K pp = K p 
and K^p = K,), and weak coupling. In this case, the 
phase diagram in the K^j K p plane is shown in Fig. |^. 



V. SUMMARY 

The main results of this paper are the interact- 
ing coulomb gas representation of the competing stripe 
and superconducting orders, Eq. (g); and the renor- 
malization group flow Eq. (17). Analysis of these 
flow equations shows that the 1+1 dimensional ver- 
sion of the theory proposed in ref. (5| supports stable 
phases corresponding to stripe order/superconducting 
order, stripe disorder/superconducting order, stripe 
order/superconducting disorder, and stripe disor- 
der/superconducting disorder. 
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APPENDIX A: DETAILS OF 
RENORMALIZATION GROUP CALCULATION 

Here we present the details for the renormalization 
of the system described by Eq. (0) and Eq. ( |s|) and 
parametrized in fourier coefficients via equations |1(J|. [TT| . 
and |lj. The details closely follow the procedure used by 
Jose et. aZ.Q. The basic idea is to introduce a small 
length scale cutoff r c , and then integrate out configura- 
tions with pairs of the same type of charge which are 
r c + dr c apart to find a new system with a longer mini- 
mum length scale. 

In our action Eq. (g), the fields N and M consist of 
point charges, e.g. N(r) — J2 a ^ a ^( r ~ Ta )- We will 
express the terms Eq. (||) involving Gn and Gm as sums 
over pairs of these charges. This results in the cancella- 
tion of all odd fourier components oi Gn and Gm, i-e. 
a n and a n are zero for odd n. To see this, note that in 
re-expressing the sum as a sum over pairs, we combine 
terms like N a NpG N (R a - Rp) + NpN a G N CR-p - R a ). 
The two Gat's above differ by reversing the relative vec- 
tor (i.e. 8 — > 8 + 7r). Since odd fourier components pick 
up a relative minus sign under this reversal, they cancel 
each other, and only the even fourier components sur- 
vive. On the other hand, since the terms with T describe 
the interaction between distinguishable vortices, we can- 
not convert them into sums over pairs, and therefore the 
fourier coefficients c n can be nonzero for both odd and 
even n. 

As in the Kosterlitz-Thouless case, Gn and Gm di- 
verge logarithmically at short length scales. In Eq. (0) 
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this translates to divergences when Ri = R2. To re- 
move this divergence, we must enforce charge neutrality 
(J2 a -^p = E Q M a = 0) and impose a small distance cut- 
off r c M. To account for the microscopic physics lost in 



this procedure, we introduce core energies Ajv and Am 
for the vortices (charges) . After this our action is written 
as 



S = N a NpG N (R a - Rp) + M a MpG M {R a - Rp) + £ M a NpiT N (R a - R^) + £ N*A N + £ M*A M 

(a, 13) (a,/3) a,f) a a 

(Al) 

where (a,/?) denotes a sum over pairs and a, (3 denotes an unrestricted sum over both a and /3. At this point we 
make the simplifying assumption that An and Am are very large, so we may restrict to N a , M a — ±1. Introducing 
the fugacities un = e _Ajv and ijm = e _Aflf , we can write the partition function as a sum over j N-dipoles and k 
M-dipoles: 

j=0fe=0 J- K - J r c r c J r c r c 

S = N a N G N (x a -x fj )+J2 M a M G M (z a ~z fj ) + J2 M a NpiT N (z a - xp) (A3) 

(a, 0) (a,/3) a,0 

The first step in the real-space renormalization procedure is to integrate over dipoles of size r c + dr c . Because 
Vn-iDm « 1 we need only consider configurations with one dipole of either type. Furthermore, mixed M-N dipoles 
are not integrated out because doing so would violate overall charge neutrality for each species. This means that the 
M and N charges are configured freely with respect to each other. Writing out these single dipole contributions, the 
partition function gains an extra factor: 

-s l\ fd 2 RN r +drC 2 -S' f d2R M r +dr *d 2 T M 2 _g, 4 , 1 

e s -» e s x I 1+ / -— - / -^Vle s - + / —f- / —^Vm* Sm +0( V 4 ) \ (A4) 

2k 

S' N = -G N (r N ) + J^N a (G N (R N + ^-x a )-G N (R N - 1 ^--x a )) + 

a 

J2 M a (iT(R N + ^L- Za )- iT(R N - T -f~ z a )) (A5) 



2j 

S' M = -GAf(rM) + ^M Q (G , M(RM + ^ 1 -z Q )-G M (RM- £ f-z Q )) + 



2 k 



N a (iT(R M + ^ - x Q ) - iT(R N - ^ - x Q )) (A6) 



2 v " 2 



where R signifies the center and r the separation of the dipole being integrated out. To obtain the above, we have 

i+i) 2 ■ 



used the fact that the combinatorial factor for j + 1 pairs of charges, one of which is a dipole (i.e. ftrrW), is equal 



to the combinatorial factor of j pairs of charges(i.e. jp) 



To proceed, rewrite equations A5 and A6 using identities like (Gn(Rn + ^f" — x a) — Gn(Rn — ^f- — x Q )) = 
?n • ^r n Gn(Rn — x a ). In Eq. (A4), these expressions enter into exponentials, which we expand. The linear terms 
in this expansion produce nothing interesting. To see this, note that there are two types of these terms: those with a 
gradient, and those without. The linear terms involving the gradient do not contribute since their integral with respect 
to r is zero. The linear terms involving no gradient (—Gn(yn) and — Gm(^m)) integrate to constants independent of 
j, k, and the N a , M a , so they merely produce an overall constant multiplying the partition function. 



Since the linear terms in the expansion of the exponentials of Eq. (A4) are uninteresting, we look at the second 
order terms. There are three types of terms here: those involving no gradients (e.g. GjyGjv); those involving one 
gradient (e.g. Gjyr ■ VG^v); and those involving two gradients (e.g. r • VGjyr • VGjy). The first type of term, with no 
gradients, integrates to a constant independent of j, k, and the N a , M Q , so we ignore it. 
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The second type of term, with only one gradient, integrates to zero. For example, one such expression is 
J2 a N a G(r N )r N ■ \7 Rn Gn(Rn — x Q ). However, in each term of the sum we can change integration variables to 
H' N = Rat — x„, producing J J2 a N a G(rN)rN • Vr^G/y (Rjy). Charge neutrality then makes the sum over N a zero. 

The third type of term, with two gradients, does contribute to renormalization. A typical term of this type is (using 
the fourier expansion of T) 
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y% M a M p r N ■ V Rn % J j!^2^2 (E 

f d 2 k 1 



c e in6 ) e ik -( z °~ R «) 



r /v ■ Vr. 



c e imS ) e ik '( z ' 3 ~ Rlv ) 



Using r/y w r c (cos(/), sin<^>) we can integrate over rjv to obtain 

fd 2 R N v-^ w f d 2 k' d 2 k . ,, , ,,_ 1 v-^ 



1 2-7T rj?<ir c 2 



inQ+imQ' ik-z^+ik'-z^ — i(k-t-k') -R/\ 



7r dr. 
2!' 



c -y 2 N Y M « M f 

a.fi 



d 2 k 1 

(2tt) 2 A; 2 



^ik-tzc-z^) 



E 



z(n+m)(9 irn-K 



it dr, 
2~ t 



c 2 



v- f d 2 ] 



d 2 k l 



(2tt) 2 /c 2 



,ik-(z a -z/3) 



E 



,<5n,/+rnC/C m (-l) 



m+1 | ^,zn# 



r 



After many calculations like this and some algebra, we 
obtain a new partition function in the form of Eq. ( A2 ) 
again, except that the interaction functions have been 
redefined to G' N , G' M , and T' via the new coefficients 

a' n = a n - nYdn,k+ m (y 2 N aka m + y 2 M {-l) m+1 c k c m )dl 

k,m 

a' n = a n -nY S n,k+ m (y 2 M a k a m +y 2 N (-l) m+1 c k c m )dl 

k : va 

c' n = c n - Trjjn^+miyjjak + y 2 M ak)c m dl. (A7) 

k,m 

Here dl = dr c /r c . 

At this point it is useful to examine more closely the 
structure of the contributions to the interaction functions 
and understand what sorts of interactions we are dealing 
with. To this end we evaluate a typical term, 



d 2 k 



(2tt) 2 k 2 
using k • r = kr cos <fr this becomes 



1 



(2n) 
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dB 



2tt 



^_ e ikr cos(9-<j)) e in9 



'-J- n {x) 



(A8) 



(A9) 



where we have introduced the Bessel function J„. Due to 
the oscillatory nature of the Bessel function the integral 
is dominated by the infrared so we employ the asymptotic 
form J n (x) ~ ^r(f )™ to obtain 



27m! ' 



incp 







| const + 


/>i (f ) 







(A10) 



where the negative sign refers to n < 0. For the case 
the infrared part converges and we are left with a 
function of the spatial angle <f> independent of r. For the 
case n = we obtain a logarithmic divergence in r. Using 
the cut-off r c the precise form for the n = contribution 
to the interaction functions is 



/ 



d 2 k e zk r 



2tt 



In 



(All) 



Thus in the interaction terms only the zero-mode 
fourier component has any mention of r c . To complete 
the renormalization group program, we must write the 
partition function in a form containing only the renor- 
malized core size r' c = (1 + dl)r c . After we do this the 
partition function is (up to an overall constant) 
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1 1 fdxi d 2 x 2l fdrzi d 2 z 2 k 

.... x 



exp(-S') x exp J --L^ln(l + dZ)i £ iV Q A^ - J-^ ln(l + ^ M a Mp - -c' ln(l + £ N a Mp 



where S 1 ' is the same as in equation A3, but with r c , Gat, 
Gm, and T replaced by r^,, G^, G' M , and T', respectively. 
The second exponential is the correction from changing 
r c to r' c in the zero modes of the interactions. 

Now, by charge neutrality N a = M a — 0. There- 
fore E a ,pM a N p = 0. Also E a ^N a N p = (J2N a ) 2 - 
Y^Nl = -2j and similarly J2 a ^/3 M a M (i = -2k, so 
that the last exponential in equation A12 becomes (1 + 

a' a' 

dl)~ 2j T^(l + dl)^ 2k T^. With these contributions the 
renormalization of y^ and yu is 



to lowest order in the y. At this point, the partition 
function is in the same form as in equation A2 except that 
r c has been replaced by r' , and the fourier coefficients and 



fuga cities have changed according to equations A7 and 
A13. This completes the renormalization group program 
and gives us the the differential renormalization group 
flow equations [l7| stated in the text. 



y' N = VN +?Mr(2- -p)dZ 

4-7T 

Vm = Vm + Vm{2 - -^-)dl 



(A13) 
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